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We study topological properties of phase transition points of topological quantum phase transitions by as¬ 
signing a topological invariant defined on a closed circle or surface surrounding the phase transition point 
in the parameter space of momentum and transition driving parameter. By applying our scheme to the Su- 
Schrieffer-Heeger model and Haldane model, we demonstrate that the topological phase transition can be well 
characterized by the defined topological invariant of the transition point, which reflects the change of topological 
invariants of topologically different phases across the phase transition point. 
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I. INTRODUCTION 

Conventional continuous quantum phase transitions (QPTs) 
are driven by pure quantum fluctuation effects due to the 
change of external parameters and generally described in 
terms of the spontaneous symmetry breaking and order pa¬ 
rameters of the ground stated. On the contrary, topological 
QPTs involve the change of ground-state topological prop¬ 
erties and accompany no symmetry breaking^.-^. Different 
topological states are classified by topological quantum num¬ 
bers, which take discrete numbers, in contrast to order param¬ 
eters used in conventional QPTs to distinguish various phases, 
which generally take continuous values. Conventionally, con¬ 
tinuous QPTs can be classified into different order QPTs by 
singularity properties of the ground-state energy at the phase 
transition point (or critical point for ^ > 2), i.e., ^th order 
QPTs are characterized by discontinuities in the ^th deriva¬ 
tive of the ground-state energy. 

As the singularity of ground-state energy plays an impor¬ 
tant role in determining universal properties around the criti¬ 
cal point of the QPT, however it can not distinguish whether 
the phase transition is a topological QPT or a conventional 
one within the Landau-Ginzburg paradigm^. Beyond the tra¬ 
ditional energy criterion, a QPT can be also witnessed by qual¬ 
itative changes of physical quantities related to the ground- 
state wavefunctions, e.g., the Berry phased, quantum fidelity 
and the fidelity susceptibility2ri2, and the quantum geomet¬ 
ric tensor—"—. Although these approaches have shed light on 
our understanding of QPTs from the geometric aspect of the 
ground-state manifold, one can not identify a QPT to be a 
topological or trivial one solely from the singularity of par¬ 
ticular physical quantities at the phase transition point unless 
additional quantities related to the topological invariant are 
calculated. An interesting question arising here is whether we 
can characterize a QPT is topological or conventional phase 
transition from the property of the phase transition point? 

To answer the question, let us recall that a topological QPT 
distinguished from a trivial one is manifested by the change 
of topological invariant, instead of symmetry breaking, across 
the transition point. While the topological invariant, e.g., the 
quantized Berry phase for one-dimensional (ID) topological 
systemsi^ or Chern number for two-dimensional (2D) quan¬ 


tum Hall systemsi^, is well defined for a gapped phase apart 
from the QPT point and characterizes the global geometrical 
property of the Bloch band, it fails to work at the gapless criti¬ 
cal point. To overcome the difficulty, in this work we propose 
an alternative definition for the topological invariant, which 
is not defined on the momentum space at the transition point, 
but via a closed detour path surrounding the critical point on 
the parameter space spanning by both the momentum and the 
transition driving parameter. By applying this idea to the ID 
and 2D topological systems, e.g., the celebrated Su-Schrieffer- 
Heeger (SSH) model and Haldane model, we demonstrate 
these topological invariants taking nontrivial quantized num¬ 
bers for topological QPTs, but some non-universal numbers 
or zero number for conventional QPTs. Our results suggest 
that we can judge topological or trivial QPTs from topologi¬ 
cal properties of the phase transition points. 

II. MODELS AND RESULTS 
A. ID topological models 

We begin our discussion with one of the simplest ID topo¬ 
logical systems, the SSH modeli^, which can be described by 
the Hamiltonian: 

H = ^[(f + S)cl .CB,i + it- + h.c., (1) 

i 

where . is the creation operator of fermion on i-th A (or 
B) sublattice. This model has two sites in a unit cell, the hop¬ 
ping amplitude in the unit cell is t S and that between two 
unit cells is t-S. For convenience, r = 1 is taken as the energy 
unit. After the Fourier transformation Csj = 
with V = A{B), the Hamiltonian can be written as 

H = (2) 

where = (c|^,c|^)and/z(A:) = hxCr^+hyCTy with cr the Pauli 
matrix acting on the vector if/k, hx = (1 + ^) + (1 - ^)cosk, 
hy = (1 - ^)sink. It is well known that this model be¬ 
longs to the BDI class according to the standard topological 
classification^^ and has two topologically distinct phases for 
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FIG. 1: (Color online) (a) The spectrum of the SSH model versus d 
under OBC. (b)-(d) The spectrum versus k under PBC, with (b) 6 = 
-0.5, (c) (5 = 0 and (d) ^ = 0.5. (e) and (g) show the winding of h(k) 
across the Brillouin zone, corresponding to (b) and (d) respectively, 
(f) shows the winding of h(6) as 6 varies a period. The arrows in 
(e)-(g) show the direction of the Hamiltonian h{k) and h{6). 


S > 0 and ^ < 0 with the phase transition point at ^ = 0, 
where the gap closes at ^ = tt. Under the open boundary con¬ 
dition (OBC), these two phases can be distinguished by the 
presence and absence of degenerate zero-mode edge statesi^, 
as shown in Fig [Ha). 

While the spectrum under periodic boundary condition 
(PBC) shows a similar structure for ^ < 0 or ^ > 0, as dis¬ 
played in Fig[TIb)-(d), the topological property of the distinct 
phase can be characterized by the Zak phase^iiS^, i.e., the 
Berry phase across the the Brillouin zone, which is defined 
as 


the Berry phase is defined as 

= i C dO((p(0)ldol(p(0)). (4) 

J-n 

After some algebras, one can obtain y^ = tt, which corre¬ 
sponds to the topological phase transition at ^ = 0. In Fig[TJf), 
we also show the winding of h(6) as 6 varies a period, giving 
rise to a winding angle of 271. 

As a comparison, we next consider a topologically trivial 
two-band model with alternating on-site potentials, described 
by the Hamiltonian: 

^ = Z ^A,i+l^BJ + h.C.] +lu(c\.CA,i - clfBj)- (5) 

i 

This model has alternating chemical potential fi and -yu for 
site A and B, and has a similar spectrum as the SSH model 
with a phase transition occurring at yu = 0. However, as yu 
breaks the inversion symmetry, the Berry phase of each band 
is no longer quantized and no degenerate edge states emerge 
under the OBC. Similarly, we can also calculate the Berry 
phase around the gap closing point at yu = 0 and k = 7i. The 
numerical result shows that yd is not a quantized invariant and 
is associated with the radius A of the integral path. On the con¬ 
trary, the yd for the SSH model is always tt regardless of the 
value of A, which suggests that yd is a topological invariant. 
This difference means that we can judge whether the QPT is 
a topological phase transition from values of the Berry phase 
yd around the critical point. 

Similar scheme can be directly applied to the other ID 
topological nontrivial systems with the ground state charac¬ 
terized by the Z-type invariant. According to the ten-fold-way 
classificationii, in one dimension the BDI (orthognal), AIII 
(unitary) and CII (symplectic) classes belong to the Z type. As 
the SSH model belongs to the BDI class, next we consider the 
fermionic Creutz Iadderi2i22, which is related to the ID AIII 
topological insulator. The Creutz ladder can be described by 
the Hamiltonian: 


y = i I dk{(f(k)\dk\cp(k)}, (3) 

d-7T 

with (f(k) the eigenstate of the occupied Bloch band. While 
the topological phase is characterized by y = tt for ^ < 0, the 
trivial phase corresponds to y = 0 for ^ > 0. The geometrical 
meaning of the Zak phase can be understood as the winding 
angle of h(k) as k varies across the Brillouin zone^^, as shown 
in Fig (He) and (g). For the topological nontrivial case with 
S = -0.5, the direction of h(k) winds an angle of Itt, whereas 
for the trivial case with 6 = 0.5 the winding angle is zero. 
However, when ^ = 0, the two bands are degenerate at k = tt 
(FigCtc)), and the Zak phase is ill-defined. 

To describe the topological property of state at the phase 
transition point, here we defined the Berry phase on a circle 
around the gap closing point in the parameter space of k and 
S. Introducing 6 as the varying angle and A as radius of the 
circle, we have k = AsinO 7i and S = A cos 6, hence the 
Hamiltonian around the circle can be represented as h(0), and 


H = 




CA,i + 


+ + h.C., 


( 6 ) 


where K and M are tunneling strengths, a is a magnetic flux, 
and ^ = 1 is set to be the energy unit. In the momentum space, 
the Hamiltonian h(k) is given by h(k) = hoi-\-hxCrx-\-hzCr^, with 
ho = 2 cos a cos k, hx = 2 cos k M, h^ik) = -2 sin a sin k. 
Generally, this model has no time-reversal symmetry due to 
the introduction of a. For the specific case with a = ±7r/2, one 
can check that the Hamiltonian fulfills the chiral symmetry 
(Tyh{k)(Ty = -h(k), which means the system with a = +7r/2 
belonging to the AIII class. This model also has a two-band 
spectrum, and the gap closes when sina = 0, M = -2 cos k 
or M = +2, cosk = +1. The phase diagram is shown in 
Fig 121 which indicates the phase within the parameter regime 
of |M| < 2 is topologically nontrivial characterized by the Zak 
phase y = +7r, whereas the phase with \M\ > 2 is a trivial 
phase with y = 0. Notice that although the Zak phase usually 
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FIG. 2: (Color online) The phase diagram of the Creutz ladder model 
with on the phase boundary being marked. There are two degen¬ 
erate points in the Brillouin zone for the system with \M\ < 2 and 
a = 7T or 0, and the summation of y^ is +2;r. 


only takes 0 or tt with a modulus 27r, the winding directions 
of h(k) are opposite to each other for y = tt and y = -tt, 
which represent different topological phases. It is interesting 
to point out that the state at a = 0 and n within \M\ < 2 is a ID 
semimetal with band touching points at ^ + arccos(-M/2). 

At the phase boundary M = +2, which separates the trivial 
and topological phases, the band gap closes at ^ = tt or 0. 
Similar to Eq.®, we can calculate the Berry phase around 
the gap closing point at (M = 2, k = n) or (M = -2, k = 
0). As displayed in FiglU we have yj = ±7t at the phase 
boundary, which corresponds to the change of the Zak phase 
y at different phases. If we fix M and take a as a driving 
parameter, we can calculate the Berry phase around the band 
touching points {a = 0,k = ± arccos(-M/2)) or (a = 7r,k = 
+ arccos(-M/2)). For a = 0, we get y^ = -tt at both gap 
closing points, and fora = 7r we get y^ = n. It is clear that the 
summation of yj for these two degenerate points gives -2n or 
271, consistent with the change of the Zak phase y at different 
sides of Of = 0 or Of = TT. 

We have indicated that the Creutz ladder system with a = 
+7r/2 belongs to the AIII class. When a deviates +7r/2, the 
chiral symmetry is broken and the system can not be clas¬ 
sified into the standard ten-fold classesii. Nevertheless, the 
system still supports topologically nontrivial phase with the 
Zak phase y = +7r as the system is protected by the inversion 
symmetry^ii2i, i.e., crxh(k)crx = h(-k). Through the above ex¬ 
amples, it is clear that our scheme works for both the standard 
topological classes, e.g., the BDI class and the AIII class, and 
the non-standard class protected by other symmetries as long 
as the topological phase can be characterized by a nontrivial 
Berry phase. 



FIG. 3: (Color online) The energy spectrum of the Haldane model 
with q = 1 and t 2 = 0.1. (a) M = 0, cr = 0; (b) M = 0, cr = 7il3\ 
(c) M = 3 V3t2 sincr, a - 7t/3\ (d) M - -3 V3t2 sincr, a - 7:13. (a), 
(c) and (d) are on the phase boundary, while (b) is in the gap opened 
region. The gap closes at both Dirac points in (a), while it closes only 
at one of the two Dirac points in (c) and (d), respectively. 

B. 2D topological models 

Next we apply a similar scheme to study the topological 
property of phase transition points of 2D systems. We be¬ 
gin our discussion with the famous Haldane model2^, which 
supports a rich phase diagram, exhibiting either topological 
or trivial phase transitions. The Haldane model is a proto¬ 
type model which may realize the anomalous quantum Hall 
effect in a 2D honeycomb lattice without any net magnetic 
flux through a unit cell of the system. The Hamiltonian of the 
Haldane model is given by: 

H ='^ tocjci + ^ ticjcj + ^ t2e‘“‘-^cjcj, ( 7 ) 

i (i,i) (iUj)) 

where the summation is defined on the 2D honeycomb lat¬ 
tice, which is composed of two sublattices labeled by A and 
B, respectively. Here to = M for site A and to = -M for site 
B, ti denotes the nearest-neighbor hopping amplitude, and t 2 
denotes the next-nearest-neighbor (NNN) hopping amplitude. 
The magnitude of the phase is set to be \aij\ = a, and the di¬ 
rection of the positive phase is clockwise, following Haldane’s 
work. This model is well known for its three topologically dif¬ 
ferent phases characterized by the Chern number C (C = +1 
or 0), which is defined as the integral of the Berry curva¬ 
ture V for each band in the Brillouin zone, with V defined 
as V X /(^(k)|V|^(k)), where ^(k) = (f(kx, ky) is the eigenstate 
of the occupied Bloch band. 

Consider the case with \t 2 /ti\ < 1/3, for which the two 
bands never overlap and only touch at the the Brillouin zone 
corner when M = +3 ^/3t2 sin a. Expanding the Hamilto¬ 
nian in the momentum space around the Dirac point K+ = 
(+^^, 0), i.e. kx = ± 5 ^ + ky = y, we get the effective 
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FIG. 4: (Color online) (a) and (b) show the phase diagram of the 
Haldane model with on the phase boundary being marked, (a) is 
for the case with defined in the space of momentum and a; (b) 
is for the case with defined in the space of momentum and M. 
(c)-(e) show the direction of /z(k) near the gap closing point in (a) 
with different ag: (c) ao = tt/S, M = 9^2/2, 0 ^ = 1 ; (d) ao = 7r/2, 
M = 3 V3t2, Q = 0 ; (e) ao = 2;r/3, M = Q = -1. 


Hamiltonian: 

h(k) = + Vs^ixcTjc + hycTy + (M + 3 ^^3t2 sin cr)(T^. (8) 

The gap closing points are at x = y = 0 and M + 3 ^^3t2 sin a = 
0. While the gap closes at both of K+ when M = 0 and a = 0 
or TT, there can be no more than one gap closing point for other 
value of a, i.e., the gap closes at K+ when M = 3 ^^3t2 sin a or 
at K_ when M = -3 ^^3t2 sin a. Different situations of band 
touching are shown in FiglJl 

The phase diagram of Haldane model is displayed in FigH] 
(a) or (b) with different phases characterized by the Chern 
number C. At the phase boundary, the gap between the up¬ 
per and lower energy bands closes, and the Chern number 
is ill-defined therein. Similar to the ID case, we can define 
a topological invariant on a closed surface surrounding the 
phase transition point in the parameter space of k and tran¬ 
sition driving parameter 6 to describe the topological property 
of the transition point. Different from the ID system, here the 
topological invariant is given by 

Q = ~^y-dS, (9) 

with S a sphere surface around the gap closing point 
(k^, ky, So). Following Berry’s work2&, the Berry curvature for 
the lower band can be expressed as: 


\ = Im 


(-|V/z(k,^)|+)x(+|V/z(k,^)|-) 
{E. - 


( 10 ) 


with l-h) (|-)) the eigenstate of the upper (lower) band. With 
some further calculations, the Berry curvature can be written 


as: 

2R^ dky dS ^ dky dS ^ dky dS ^ 

dh^dhy 

dky dS ^ dky dS ^ dky dS ^ ’ 

J_ . dh^dh^ ^^ 

2R^ dS dkx ^ dS dkx ^ dS dkx ^ 

dhz dhy ^ 

dz dkx ^ dz dkx ^ dz dkx ^ ’ 

1 dhx dhy dhy dh., dh. dhx, 

2R^ dkx dky ^ dkx dky ^ dkx dky ^ 

^ 

dkx dky ^ dkx dky ^ dkx dky ^ 

Here the Hamiltonian h(k) = hxCTx + hyCTy + and R = 
hih^ hj. For convenience, we choose a spherical sur¬ 
face surrounding the gap closing point with a radius of A, i.e., 
kx = kl X, ky = ky y, S = So z with x = A sin ^ cos 0, 
y = A sin ^ sin 0 and z = A cos 6, where 6 is the polar angle and 
0 is the azimuthal angle of the spherical surface. The integral 
then becomes: 


= - 


y^ = - 


y" = - 



(sin 6 cos 0, sin 6 sin 0, cos 0)A^ sin 6d6d(p. 


If all the hx, hy and h^ are linear of (v,y,z) around a gap 
closing point, we can rotate and stretch the axes to reach 
hx = ±x, hy = +y, hz = ±z, and the Chern number Cd = ± 1. 

As the phase boundary is given by \M\ = 3 ^/3\t2 sinal, ei¬ 
ther the parameter a or M can be chosen as the phase transi¬ 
tion driving parameter. Similar to the ID cases, we can judge 
whether the phase transition is topological or not by examin¬ 
ing Cd defined around the gap closed points. First we choose 
a as the third parameter besides the momentum k by keeping 
M fixed. We take the case with the transition point located at 
K+ as an example, and the case at K_ can be analyzed sim¬ 
ilarly. Defining a = ao + z with M - 3 ^^3t2 sin ao = 0 and 
expanding the Hamiltonian near ao, we have 


h(k) = - ^|3t\xo■x + hyo-y - 3 ^f3t2 cos aozcz (11) 


when \M\ < 3 V3|^2l, and the Chern number around K+ is 
Cj(a) = sgn(cosao). When \M\ = 3 ^^3\t2l i.e., the case with 
ao = ±7112 marked by the “cross” in FigUa), however, the 
expansion becomes 


h(k) = - ^[3tlXC^x + ^lycr^ + 3 ^^3t2Z^c^z, (12) 


and the integral results in Cj(a) = 0 as the integrand is an odd 
function in the interval. In Fig Ha), we show the value of Cd 
around the phase boundary: except Cj = 0 at the four points 
marked by “crosses”, it takes either I or -I. We can see that 
Cj(a) shows the change of the band Chern number C across 
the transition point by varying a. For the case of Cd = 0, the 
change of C is zero when varying a, which indicates a topo¬ 
logically trivial phase transition. On the other hand, the case 
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FIG. 5: (Color online) (a)-(c) The energy spectrum of the QWZ 
model with p - 1, (a) yu = -2, (b) yu = 0 and (c) jd - 2, respec¬ 
tively. (d)-(f) display the direction of the Hamiltonian h in the in the 
x-y plane, corresponding to the gap closing points in (a)-(c) respec¬ 
tively. In (e) we only display the direction of Hamiltonian around 
one of gap closing points at {k^ - n.ky - 0) shown in (b), while the 
other one at {k^ - 0,ky - n) has exactly opposite direction in both 
the V and y direction. 


of Cd = ±l corresponds to a topological phase transition from 
the trivial (topological) phase to topological (trivial) phase. 
Particularly, when M = 0 and a = 0 (or n), there are two gap 
closing points in the Brillouin zone, and the change of C is +2 
when varying a, which corresponds to the summation of Cd 
around these two points. 

If we choose M as the third parameter besides the mo¬ 
mentum k by keeping a fixed and define M = Mq z with 
Mq + 3 V3 sin cr = 0, the expansion of the Hamiltonian be¬ 
comes: 

h(k) = + ^^3tlxcrx + hycTy zcr^, (13) 

and the Chern number around the gap closing point of K+(K_) 
is always Cd(M) = -1(1), as shown in FigUb). The value of 
Cd(M) also indicates the change of C when varying M. When 
M = 0 and a = f) and tt, the change of C is zero, which also 
matches the summation of Cd(M). 

To visualize the Chern number Cd, we show the direction of 
h(k) around the gap closing point with different ao in FigHJc)- 
(e). One can see that all the arrows point to z = 0 in (c), and 
to the opposite direction in (e), but the z-component of h(k) in 
(d) is always positive. The value of Cd is associated with the 
direction of h(k). In FiglUc), the direction of h(k) is always 
toward x = z = 0 and away from y = 0, which corresponds 
Cd = 1. In FiglUe), the direction of h(k) is opposite to the 
one in FigHtc) only in the z direction, corresponding to Cd = 
- 1. However, in FiglUd), the z-component of h(k) follows the 
same direction, corresponding to Cd = 0. 

Our theory can be applied to study other 2D Z-type topo¬ 
logical systems with the topological invariant characterized by 
the Chern number. To give an additional example, we next in¬ 
vestigate a lattice model on a square lattice described by the 
Hamiltonian H = h'cr with cr = (ctx, CTy, cr^) the Pauli matrices 
and h = (sink^, sinky,p[ju - coskx - cosky]), which is known 
as the Qi-Wu-Zhang (QWZ) model^i. Depending on the value 
of Id, this model has three different topological phases charac¬ 


terised by the band Chern number C. Setting p = 1, the band 
Chern number C = 0 while \jd\ > 2, and C = -sgn(jd) when 
\jd\ < |2|. There are three phase transition points for this model: 
Id = -2, yu = 0 and n = 2. As shown in Figl3a)-(c), when 
Id = ±2, there’s only one gap closing point (kx = ky = n for 
Id = -2 and kx = ky = 0 for id = 2) in the Brillouin zone, but 
when Id = 0, there are two gap closing points (kx = n, ky = 0 
and kx = 0, ky = n). By choosing id as the third parameter 
beside the momentum, the expansion of the Hamiltonian near 
the gap closing points results in 

h(k) = -xcTx - ycTy -f zcTz, ii = -2; 

h(k) = ±xcrx + ycr3; -f zcTz, n = 0; 

h(k) = xcTx -f ycr^; -r zcr^, id = 2. 

In Figl3d)-(f) we show the direction of h in the x - y plane, 
as the z direction of h(k) is always away from z = 0. The 

result shows that Cd = 1 for yu = ±2, and Cj = -1 for both 

gap closing points when id = 0, hence the summation of Cd at 
each phase transition point also corresponds to the the change 
of the band Chern number C. 


III. SUMMARY AND OUTLOOK 

In summary, we proposed a scheme to study the topolog¬ 
ical properties of phase transition point for various topolog¬ 
ical QPTs by introducing a topological invariant defined on 
a closed curve in the parameter space surrounding the tran¬ 
sition point. By studying several typical topological models, 
we demonstrated that the topological or trivial phase transition 
can be distinguished by the introduced topological invariant 
around the transition point, which takes nontrivial quantized 
numbers for topological QPTs but non-universal numbers or 
zero number for conventional QPTs. Our theory provides a 
way to discriminate topological or trivial QPTs by directly 
studying the properties of the phase transition point. 

In order to give a proper definition to the topological in¬ 
variant of a phase transition point, we have introduced the 
phase transition driving parameter as an additional dimension 
parameter besides the momentum space, and thus raise the 
effective parameter dimension by one. Such an idea seems 
to share some similarities with the construction in the anal¬ 
ysis of topological phases by dimension extension through 
introducing an additional dynamical parameteri22rii. While 
the topological invariant is defined in the enlarged (D 1)- 
dimensional parameter space for D-dimensional systems in 
previous works^Szi^, in the present work the topological in¬ 
variant is defined on a close curve of (D l)-dimensional 
space surrounding the phase transition point and thus is ef¬ 
fectively defined on a D-dimensional curve. This is the main 
difference from the previous construction via dimension ex¬ 
tension. While our scheme can be applied to Z-type topolog¬ 
ical systems, it can not be directly generalized to deal with 
the Z 2 -type systems as the corresponding topological invari¬ 
ant is described by the Z 2 number—. It would be interesting to 
study how to characterize the topological phase transition for 
^ 2 -type systems by analysing the topological property of the 
phase transition point in the future work. 
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